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GRB 990123 revisited: Further Evidence for a Reverse Shock 
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ABSTRACT 



Recently we have presented a new theoretical analysis of the reverse shock 



O 

o 

C\| ■ emission. We use this analysis here to revisit GRB 990123. We find new and 

I>y compelling evidences that the optical flash and the radio flare of GRB 990123 

resulted from a reverse shock. This suggests that a significant fraction of the en- 
ergy of the relativistic ejecta must have been carried by baryons. It also suggests 
that the external medium is an ISM and that in this burst the reverse shock 
emission dominated at early time over other possible processes. We use the early 
optical emission to constrain the physical parameters of the original ejecta and 
the microscopic parameters in the emitting reverse shocked region. 
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Introduction 



The 9th magnitude optical flash of GRB 990123 was one of the most exciting discoveries 
associated with Gamma- Ray Bursts (GRBs). This was followed by a detection of a strong 
unusual radio flare from the same source (Kulkarni et al. 1999, Galama et al. 1999). One of 
the most impressive facts concerning the strong optical flash of GRB 990123 was that such 
a strong prompt optical emission was predicted, just a few months earlier to arise from the 
^ ■ reverse shock during the early afterglow (Sari, & Piran 1999a). Later this interpretation was 
questioned and other models for prompt optical emission were proposed (e.g. Beloborodov, 
2002). We return to this issue here. 

In the popular fireball model (Piran 1999) for Gamma- Ray Bursts (GRBs), the afterglow 
results from a blast-wave that propagates into the circum-burst medium. This blast-wave is 
originated by the energy dissipation of the relativistic ejecta to the circum-burst medium. 
In this model the early afterglow is produced during this dissipation, and thus can be used 
to study the nature of the ejecta. In case that the ejecta is Baryonic, a reverse shock is 
produced during this dissipation. This reverse shock is predicted to produce strong optical 
flash and radio flare (Sari & Piran 1999b, c). 
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Recently we have presented new theoretical analysis of the reverse shock emission (Nakar 
& Piran 2004; hereafter NP04). This analysis includes a new test of reverse shock emission 
and diagnostic tools for the physical parameters of the original ejecta. We use these new 
tools, here, to re-investigate the early afterglow of GRB 990123. Note that numerous authors 
analyzed the early afterglow of GRB 990123 (Sari & Piran 1999b; Meszaros & Rees 1999; 
Kulkarni et al. 1999; Galama et al. 1999; Wang, Dai & Lu 2000; Panaitescu & Kumar 2001; 
Fan et al. 2002; Soderberg & Ramirez- Ruiz 2002). We don't attempt to review here these 
results and compare them as the scope of such a review is too long for a short letter. 

The reverse shock optical emission has a typical and characteristic t~ 2 decay that was 
observed in the light curve of GRB 990123 (Fenimore et al. 1999; Sari & Piran 1999b, c). 
We find additional conclusive evidence that the optical flash and the radio flare of GRB 
990123 indeed resulted from the emission of a reverse shock produced by an interaction of 
the ejecta with and ISM. These results suggest first that at least in this burst the reverse 
shock emission, and not another mechanism such as a pair enriched forward shock (due 
to interaction of the prompt 7-rays with the circum burst medium; Thompson & Madau 
2000, Beloborodov, 2002), dominated the early afterglow. As the reverse shock requires, a 
baryonic component in the ejecta it also implies that a significant fraction of the energy of 
the relativistic ejecta was carried by baryons. 

The sparse observations before the peak of the optical flash (only one measurement 
before the peak at each wavelength) poses difficulties in analyzing the properties of the 
ejecta, since this analysis depends on the rising power-law index of the optical light curve 1 . 
Still, we are able to determine that the reverse shock was mildly relativistic and that the 
initial Lorentz factor is most likely smaller than ~ 600(-E 54 /n) 1 / 8 , where £"54 is the (isotropic 
equivalent) kinetic energy in the ejecta after the prompt emission phase in units of 10 54 ergs 
and n is the protons density of the circum-burst interstellar medium in c.g.s. The width 
of the ejecta is found to be, as expected, roughly equal to the burst duration (multiplied 
by the light speed). Because of the sparse data available we obtain only upper limits to 
the initial Lorentz factor and the shell's width and only lower limits to the microscopic 
equipartition parameters and the external density. We expect that in the future detailed 
optical light curves, provided by Swift and rapid follow up observations, would enable an 



1 Since there is only a single observation before the peak, the rising phase may also be very irregular 
with no power-law behavior. In this case the analysis of the ejecta properties presented in NP04 does not 
apply as it imply that the relativistic ejecta hydrodynamical profile is highly irregular. However, having no 
detailed observations of this phase we use here the simplest light curve that fits the observations - a broken 
power-law. Note that the tests, which show that the emission is a reverse shock emission, are independent 
of the shape of the rising light curve, and can be carried also if it is irregular. 
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accurate determination of the initial physical parameters of the relativistic wind by using 
the method we apply here. 

For completeness we begin (§2) with a short summary of the theoretical model of the 
reverse shock emission presented in NP04. We summarize, in §3, the early afterglow obser- 
vations of GRB 990123. The analysis of the optical light curve and the resulting constraints 
on the physical parameters of the relativistic wind is presented in §4. In §5 we test whether 
the observations fit with a reverse shock emission. The results, their implications to GRB 
modelling and future prospects are summarized in §6. 



2. The reverse shock emission - Theory 

We begin with a brief summary of the main results of NP04. These results give the 
theoretical predictions of the optical and the radio emission of a baryonic reverse shock in 
an ISM environment. NP04 suggest to model the early optical light curve (similarity to 
late afterglow parametrization; Beuermann et al. 1999) as a broken power-law with five 
parameters: 

where a± > and a 2 < are the asymptotic power-law indices, F^ opt (t ) = Fq is the peak 
flux and the parameter s determines the sharpness of the peak (a large s corresponds to a 
sharp peak). The physical parameters which determine the emission of the reverse shock are: 
(1) The (isotropic equivalent) energy E, the width A and the initial Lorentz factor T of the 
ejected wind. (2) The circum burst density, n, that we consider as a constant interstellar 
medium (ISM). (3) The microscopic parameters in the emitting region, namely the energy 
equipartition parameters e e and e# that describe the ratio of the electrons and Magnetic 
field energy to the total internal energy, and p, the power-law index of the electrons' energy 
distribution. 

The strength of the reverse shock is determined by the dimensionless parameter £ (Sari 
& Piran 1995): 

I 1 / 2 

where / = (3E / (AimrripC 2 )) 1 ^ (m p c 2 is the proton rest mass energy). According to this 
definition: 

r = i88r 3/4 Ar 2 3/8 (iWn) 1/8 . (3) 
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The peak of the optical emission is observed at 

to = ^(l + 0.7e 5/2 )(l + ^), (4) 

when the reverse shock finishes crossing the ejecta. c is the speed of light and z is the redshift. 
The structure of the optical and the radio reverse shock light curves depend on the relative 
values of the three reverse shock break frequencies at t : u r a1 the self-absorbtion frequency, 
v r m , the synchrotron frequency and u r c the cooling frequency. In principal numerous patterns 
are possible. However, over a wide range of the parameter space (and as we show below, 
in particularly for GRB 990123) these frequencies satisfy u ra( n < ^m(^o) < ^a(^o) < v ovt < 
u^ito), where v r adio and u opt are the observed radio and optical frequencies respectively. We 
consider this frequencies sequence as the generic case. In this case the peak optical flux is: 

r = 0.1mJy(l + ,)-^1.5 2 - 5 - (^f) P " elZ\ef_ 2 n^Elth-jD^A^ (5) 

where throughout the paper we denote by Q x the value of the quantity Q in units of 1Q X 
(c.g.s), and D is the proper distance. The function A r F0 (^) is approximated in the range of 
0.1 < f < 2.5 by: 

^,o(0^180e a65 (6-10- 4 r 2 - 6 )^. (6) 
A r F0 (£) values out of this range can be found in NP04. 

The flux before to depends strongly on £, and thus the rising slope, a\ can be used to 
determine £ using (in the generic case): 

1.2(0.5 + |(£-0.07£ 2 )). (7) 

The reserve shock decay index is, however, a constant: 

d2 ~ -2, (8) 

and its main use is to identify the reverse shock emission. Finally, the sharpness parameter, 
s, depends mainly on the profile of the ejecta. A sharp break with s > 3 is expected when the 
ejecta is rather homogenous, while a gradual break with a lower value of s ~ 1 is expected 
if the spread of the initial Lorentz factor within the shell is large. 

In contrast to the optical emission, the radio continues to rise at t > t and it peaks at 
a later time, t*, when v Ta dio = v r a . Over a wide range of £ values (when the shock is not ultra 
relativistic) v r adio < ^(^o) < ^(^o) ~ 10 12_13 Hz. In this case the radio emission is expected 
to rise first as ~ t 0,5 until Vradio = v r m and then as ~ t L25 until t*. At t > t* it is expected 
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to decay, similarly to the optical emission, as ~ t 2 . This behavior provides a second and 
independent test for the reverse shock emission: 

p— 1 | ^ g p — 1 

y(t) 2 =f^ L ) 2 ~ 1000 > ( ) 

-TO \ r 0/ \ v radio J 

where this value can be larger or smaller by a factor of ~ 3 (for a given p) due to uncer- 
tainty in the hydrodynamics (see Kobayashi & Sari 2000). It provides also an independent 
measurement of Va(t ): 

to 

Detailed radio observations at t < t* that identify the break during the rising phase, when 
u m = u radio, would enable determination of ^(t ) as well. 



3. Observations 

At z — 1.6 GRB 990123 is one of the brightest 7-rays bursts observed so far (Briggs et 
al. 1999). The main event of this burst is composed of two very energetic and hard pulses 
which last together ~ 25sec. These two pulses contain almost all the emitted energy during 
the burst. A much softer and less energetic emission (< 120keV) is observed ~ 20sec before 
and up to ~ 50sec after this event. The total isotropic equivalent energy emitted in 7-rays 
during the burst is ~ 10 54 ergs. The timing of the main 7-rays event differ between different 
detectors. While BATSE triggered on the first soft emission, ~ 20sec before the main event, 
BeppoSAX triggered 18sec later than BATSE on the beginning of the main event. 

GRB 990123 is also the only burst where optical emission was observed during the 
prompt emission (Akerlof et al. 1999). The first detection of 11.7mag is recorded ~ 7sec 
after BeppoSAX trigger. The optical emission peaks at the second snapshot that took place 
~ 32sec after BeppoSAX trigger with 8.9mag (~ 0.8Jy)! After the peak, a fast decay is 
observed with a power law slope of « —2 which becomes shallower at later times. 

The early radio observations of GRB 990123 show also an unusual flare (Kulkarni et 
al. 1999). The first observation after 0.25days shows a 8.46GHz intensity of 62 ± 32/iJy. 
The flux in this band rises to a peak flux of 260 ± 32/iJy after 1.25days and then decline 
rapidly. Observations at other wavelengths from this epoch (~ lday) (Galama et al. 1999) 
show a flux of 118 ± 40;uJy at 4.88GHz and upper limits of several hundreds /zJy in 15, 86 
and 222GHz. 
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4. Constraining the burst parameters 

Already in 1999 Fenimore et al. (1999) pointed out that the initial decay index of the 
optical flash a 2 ~ —2 agrees with the predictions of a reverse shock emission (Sari & Piran 
1999a,b,c). At later times the decay becomes more moderate as expected from the forward 
shock contribution. Motivated by this observation, we first use the optical observations to 
constrain the physical parameters of the burst assuming that the flash arises from a reverse 
shock emission. Later, we carry the additional tests of reverse shock emission, described in 
NP04. 

When modelling the early afterglow, it is important to find the time when the main part 
of the relativistic ejecta is emitted from the source. It is not necessarily the trigger time. 
In GRB 990123 the radiated energy, and therefore most likely the energy ejected from the 
source, is clearly dominated by the main two pulses. Kobayashi, Sari & Piran 1997 (see also 
Nakar & Piran 2002), have shown that in internal shocks the observed time of the 7-rays 
pulses reflects the emission time of the wind from the source. Therefore we estimate that 
the main part of the relativistic wind is beginning to be ejected by the source at the time 
that the first dominant pulse starts rising - the BeppoSax triggering time (18sec after the 
BATSE trigger). This is the point where we set the observer clock to t — 0. According 
to this setting the duration of the burst is ~ 25sec and the peak of the optical emission is 
observed at t — 32sec. 

The observations of the optical flash of GRB 990123 (Akerlof et al. 1999) include only 
one data point before the peak. Since there are two observations after the peak which show 
similar decay (t = 57sec and t = 142sec), a 2 ~ —2, a broken power law fit is possible only if 
the peak is between the first and the second observations. Thus, we have only lower limits 
to F > 0.8Jy, cci > 2, an upper limit to to < 32sec and s can take any value (see Fig. 1). 
As we have only a lower limit on a±, Eq. (8) yields only a lower limit on £ (we use p = 2.3 
thorough out.): 

£>1- (11) 

Eqs. (4 & 3) result in the upper limits: 

A < 2.5 • 10 n cm, (12) 

and 

F \ V8 



r G < 600 . (13) 

It is reassuring to find that the width of the shell obtained for £ ~ 1 divided by the speed of 
light is similar to the duration of the GRB (~ lOsec in the burst's frame). 
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Next we use Eq. (5) to constrain the microscopic equipartition parameters, e e and e#, 
and the external density, n. We consider the total energy emitted in 7-rays, ~ 10 54 ergs, as a 
reasonable estimate of the remaining energy in the ejecta, E. Therefore, we do not attempt 
to estimate it from the optical flash observations. Instead we express the dependance of the 
resulting constraints on the value of E. Taking £ > 1 and requiring F > 0.8Jy leads, using 
Eq. (5), to: 

S 15^ <I+ " 8 » (14) 

Another constraint follows the requirement that v opt < v r c {to) (otherwise the optical light 
curve would not show the generic behavior of a<i ~ 2) 2 : 

n4\ < 20^4 1/2 . (15) 

These two constraints (together with the trivial one, e e + e# < 1) result in the following 
limits: 

e e > 0.1£ 5 - 4 - 75 , (16) 
5 • l^E^n-^ <e B < O.li^V 2 / 3 , (17) 
n > 5 • lQ- 3 E^ A cnT 3 (18) 



5. Was the optical flash a result of a reverse shock emission? 

We have already mentioned that the optical flash of GRB 990123 passes the first test 
of a reverse shock emission: ~ —2 (Fenimore et al. 1999, Sari & Piran 1999a,b,c). NP04 
have shown that in the generic case [y ra dio < v r m < u r a < u opt < v c at to) the reverse shock 
emission results also in a tight relation between the radio and the optical emission, Eq. (9). 
For the observed frequencies of GRB 990123 (u opt = 5 • 10 14 Hz and u radio = 8.6GHz) and 
p = 2.3: 

p-i 

(^EL\ 2 = 1300, (19) 

V V r adio J 

The observations of GRB 990123 show a remarkable agreement with this prediction (taking 
the times and the fluxes of the flash and the flare as the times and the fluxes of the peak 
observation, t = 32sec, F = 0.8Jy, t* = 1.25day and F* = 260/iJy): 

W /i \ £ 2 i+1 - 3 



- 800 - 4000, (20) 

Fq \to / 



2 The requirement that v radio < ^(*o) < v opt is satisfied for any reasonable value of the parameters and 
therefore cannot be used as a constraint. The values of ^(^o) an d v c(to) can ^ e found in NP04. 
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where the result include the uncertainty in the hydrodynamics. NP04 show also that if the 
above test is passed, then the self absorption frequency at t , u^(tO), can be estimated from 
the the ratio between t* and t : 

<(to) ~ ^radio = 3 ' 10 13 HZ. (21) 

This frequency is determined by the physical conditions at to- As a consistency check we 
compare this result to the possible values of v r a {to) obtained by the limits on the physical 
parameters determined from the analysis of the optical emission. We consider the case of 
£ = 1 and p=2.3 for which: 



/. I (i „ i £±2 _ N p+6 _3p+10 



p+4 



v r a (t ) = 6 • 10 u Hz [(1 + z)-—e e p -[e B ^ 2 (nE 54 )—t- 2 8 j . (22) 

Within the parameters space of e e ,e B and n which satisfy Eqs. (14 & 15) and taking E 54 = 1 
we obtain 1.5 • 10 13 Hz < z^(to) < 3.5 • lO^^^Hzl Although the dependance of this range 
on the value of E is not trivial, it is weak and the value v r a {to) = 3 • 10 13 Hz is consistent 
with the allowed range for any reasonable value of E, 0.5 < -E54 < 10. This constraint on 
the value of ^(to) is independent of the value found in Eq. (21). Once more the agreement 
between the two estimates is remarkable. 

So far we focused on the radio observations at 8.46Ghz. The lower radio fluxes both 
above and below 8.46GHz, at t*, are consistent with reverse shock emission, which predicts 
that the peak frequency at t* is the observed v Tadio (8.46Ghz in this case). Finally we note 
that although radio observations from t < can theoretically be used to farther constrain the 
conditions at to, the large errors of the single radio observation from this epoch is insufficient 
to do so. 



6. Conclusions 

We have analyzed the early afterglow of GRB 990123 according to the theoretical pre- 
dictions of the reverse shock emission presented in NP04. We find that GRB 990123 shows 
the clear signature of a reverse shock emission in two independent and robust tests. Apart 
from these two tests we carry a consistency check between the radio emission and our analysis 
of the optical emission which is passed successfully. This consistency check gives us further 
confidence both in the analysis and in the fact that the optical flash and the radio flare of 
GRB 990123 are a generic case of a reverse shock emission. 

The consistency with the reverse shock model further suggests, but does not proof of 
course, that the optical flash of GRB 990123 did not arise from a pair loaded forward shock 
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(Beloborodov, 2002). In turn this indicates that the ejecta of GRB 990123 had (at least at 
this stage) a significant baryonic component, as otherwise, for example in a pure Poynting 
flux flow, a reverse shock is not expected. 

Our analysis of the optical emission shows that the reverse shock was mildly relativistic 
and that the initial Lorentz factor is less than 600(-E 54 /n) 1 / 8 . We find that the width of 
the initial shell is similar to the duration of the burst or smaller. We constrain also the 
equipartition parameters and find that e e > O.IE^ while 5 ■ 10~ 3 E^" 6 n~ 2 ^ 3 < e# < 
0.1E^ 3 n~ 2 / 3 . The external density is limited by n > 5 • lO^E^^cm -3 . 

The determination of the initial conditions of GRB 990123 and the constraints of the 
microscopic parameters and the external density is very limited due to the sparse optical 
and radio observations. Swifts and its follow up observations would hopefully provide during 
the next few years detailed light curves. These would enable a more detailed analysis and a 
much better determination of the parameters of the burst's relativistic outflow. 
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Fig. 1. — Possible light curves (with the value of s beside each light curve) fitted according 
to Eq. (1) to the optical observations (full dots) of GRB 990123. 



